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Abstract 

, It has been recognized recently that fractional calculus is useful for handling scaling struc- 

tures and processes. We begin this survey by pointing out the relevance of the subject to 
physical situations. Then the essential definitions and formulae from fractional calculus are 
summarized and their immediate use in the study of scaling in physical systems is given. This 
is followed by a brief summary of classical results. The main theme of the review rests on the 
notion of local fractional derivatives. There is a direct connection between local fractional dif- 
ferentiability properties and the dimensions/ local Holder exponents of nowhere differentiable 
functions. It is argued that local fractional derivatives provide a powerful tool to analyse the 
pointwise behaviour of irregular signals and functions. 
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1 Introduction 

In the 1880's, contrary to the notion existing till then that a continuous function must be differen- 
tiable at least at some point, examples were constructed |IJ to demonstrate that it is possible for 
continuous functions to be no where differentiable. Weierstrass's construction provided one such 
'"q \ early example. For about three decades after the construction of such functions, they were still 
considered to be rather pathological cases without any practical relevance to the physical world. 
Perrin was the first to point out their application in a real physical situation, viz., the problem of 
Brownian motion. He also prophesied |§ the possibility of a simpler way of handling these functions 
than differentiable ones. 

Recent developments in nonlinear and nonequilibrium phenomena suggest that such irregular 
functions occur much more naturally and abundantly in formulations of physical theories. The work 
on Brownian motion (see showed that the graphs of projections of Brownian paths are nowhere 
differentiable and have dimension 3/2. A generalization of Brownian motion, called fractional 
Brownian motion, [|, is known to give rise to graphs having dimension between 1 and 2. It 
was also observed that typical Feynmann paths || §, like the Brownian paths, are continuous but 
nowhere differentiable. In fluid systems, passive scalars advected by a turbulent fluid have been 
shown J7|, § to have isoscalar surfaces which are highly irregular, in the limit of diffusion constant 
going to zero. Also there exist situations in which one has to solve a partial differential equation 
subject to fractal boundary conditions, such as Laplace equation near a fractal conducting surface. 
For instance, as noted in reference H irregular boundaries may appear to be non-differentiable 
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everywhere upto a certain spatial resolution, and/or may exhibit convolutions over many length 
scales. Nowhere differentiable functions can be used to model these boundaries. 

In the study of dynamical systems theory, attractors of some systems are found to be continuous 
but nowhere differentiable ||10|| . We will consider one specific example [10, |Tl| of the dynamical 
system which gives rise to such an attracter. Consider the following map. 



Vn+1 
z n+l 



2x n + y n modi, 
x n + y n mod 1, 
\z n + cos(27rx n ) 



(1) 



where x and y are taken mod 1 and z can be any real number. In order to keep z bounded, A is 
chosen between and 1. The equations for x and y are independent of z. The x — y dynamics 
are chaotic and are unaffected by the value of z. It was shown in [10] that if A > 2/(3 + v^>), the 
attracter of this map is nowhere differentiable torus. One more example in 2-dim phase space is 
given in section ||. 

All these irregular functions are best characterized locally by a Holder exponent. We will use 
following general definition of the Holder exponent h(y) which has been used by various authors ]12|, 
|13|1 recently. The exponent h(y) of a function / at y is given by h such that there exists a polynomial 
P n (x) of order n, where n is the largest integer smaller than h, that satisfies 



\f(x)-P n (x-y)\ = 0(\x-y\ h ), 



(2) 



for x in the neighbourhood of y. This definition serves to classify the behavior of the function at y. 

Mult if r act al measures have been object of many investigations [II, 15, [16|, [T7|, [18|, 19, 20, ^1 . 
This formalism has met with many applications. Its importance also stems from the fact that 
such measures are natural measures to be used in the analysis of many phenomenon [^2], p3"| . The 
first mathematical rigorous results concerning multifractals were given in p4j| . It may however 
happen that the object one wants to understand is a function (e.g. a fractal or multifractal signal) 
rather than a set or a measure. For instance one would like to characterize the velocity field of 
fully developed turbulence. Studies on fluid turbulence have shown existence of multifractality in 
velocity fields of a turbulent fluid at low viscosity [|14] . 



In the case of functions having same Holder exponent h, with < h < 1, at every point ( 
Weierstrass function (section |5]) is such a example) it is well known [|j] that the box dimension 
of its graph is 2 — h. On the other hand, there are functions f(x) which do not have constant 
exponent h at every point but have a range of such exponents. The set of points having same 
exponents {x \ h(x) = h} for such a function may even constitute a fractal set. In such situations 
corresponding functions f(x) are multifractal. 

In the light of these findings there has been a renewed interest, chiefly among mathematicians, 
in studying pointwise behaviour of multifractal functions. New tools have been developed to study 
local behaviour of functions, which will be of practical use. Notable among them is the use of 
Wavelet transforms 
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which have been used with some success to this effect. This 
transform allows one to study the behaviour of functions at different scales and hence it is popularly 
known as a mathematical microscope. Some well-known functions have been reanalyzed using these 
techniques |29|, |30], [H], |32[ and shown to possess a spectrum of Holder exponents, thereby establishing 
their multifractality. 
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In this article we review a more direct method, which uses fractional calculus formalism, to 



characterize local behaviour of functions. The fractional calculus formalism |34|, |35[ allows taking 
derivatives of fractional order. However, as pointed out in section |2| below, there is a multiplicity of 
definitions of fractional derivatives. Hence it is important to recognize the appropriate one which is 
suitable for extracting the scaling behaviour. Recently it was realized j36| that the concept of local 
fractional derivative, introduced in section |] below in a more general setting, is ideally suited for 
this purpose and that there is a direct quantitative connection between the lack of differentiability 
of a fractal function and the dimension of its graph. Further it provides an alternative way of 
characterizing the local Holder exponent. In this formalism, unlike in wavelet transform, Holder 
exponent shows up naturally since the magnitude of local fractional derivative (LFD) jumps when its 
order crosses a critical order. The main theme of this review is centered around this developement. 

The organization of the paper is as follows. In section II we present a brief review of fractional 
calculus formalism with its recent applications. The next section discusses fractional differentiability 
using Weyl's definition and certain classical results. In the fourth section we present our definition 
of LFD in more general setting and its connection to generalized Taylor series. In section V we apply 
this definition to a specific example, viz, Weierstrass' nowhere differentiable function. It is shown 
that this function, at every point, is locally fractionally differentiable for all orders below 2 — s, but 
not so for orders between 2 — s and 1, where s, 1 < s < 2, is the box dimension of the graph of 
the function. In the next section a general result has been proved, showing relation between local 
fractional differentiability and local Holder exponent. In section V we demonstrate use of LFD in 
unmasking isolated singularities and in the study of pointwise behavior of multifractal functions. 
Section VI concludes the article, pointing out a few possible consequencs of our results. 



2 Fractional calculus and scaling phenomenon 

Fractional calculus [^, [33], |35| is a study which deals with generalization of differentiation and 
integration to fractional orders. There are number of ways (not neccesarily equivalent) of defining 
fractional derivatives and integrations. We list here some of them, which will be used in this work. 
We begin by recalling the Riemann-Liouville definition of fractional integral of a real function, which 
is given by |33], (35f 



<Pf(x) _ 1 
[d(x — a)] q T(~q)Ja (x — y) 



f(y) 



—-rdy for q < 



0. 



where the lower limit a is some real number and of the fractional derivative 



d q f(x) 



d" 



f(y) 



[d(x — a)] q T(n — q)dx n Ja [x 
Fractional derivative of a simple function f(x) 

d g x p r(p + i) 



y) 



q—n+1 



dy for n — 1 < q < n. 



p > — 1 is given by |J3], |3"5 | 
-x p ~ q for p > —1. 



dx q Y[p — q + l) 

For future use we note the following properties of fractional derivatives. For < q < 1, 

d q sin(x) d q ~ 1 cos(x) 



dx q 



dxi- 1 



(3) 



(4) 



(5) 



(6) 
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and 



d q cos(x) 



x 



d q 1 sin(x) 



dxi T(l-q) dxi- 1 
Further the fractional derivative has the property (see ref 



viz, 



d q fWx) m d q f{(3x) 



(7) 



(8) 



dx q d(/3x)i 

which makes it suitable for the study of scaling. Note the nonlocal character of the fractional 
derivative and integral in the definitions equations (If) and (0). Also it is clear from equation (j|) 
that unlike in the case of integer derivatives the fractional derivative of a constant is not zero in 
general. 

Since the definition of fractional derivative contains a lower limit 'a', it is natural that different 
definitions of a fractional differentiability would arise depending on the lower limit one chooses. 
Weyl defined fractional derivatives, with the arbitrary limit a in (|) going to — oo, as follows. 

1 d n r f(y) 



D q f{x) 



-dy for 



n 



1 < q < n. 



(9) 



r(n - q) dx n J-oo (x - y y- n + i 

The merit of this definition is that the fractional derivative of a periodic function according to this 
definition is a periodic function, which is why it is suitable in harmonic analysis. 



2.1 Applications to the study scaling in physical systems 

Recent work shows that the fractional calculus formalism is useful in dealing with scaling phenom- 
ena. The purpose of this subsection is to point out a few relevant parallel developments even though 
they are not the part of the central theme of this paper. 

Hilfer |]3~T , |38| , [39] has generalized the Eherenfest's classification of phase transition using the 
fractional derivative. Some recent papers [[K], 41], ^2], |43| indicate a connection between fractional 
calculus and fractal structure |3], ||] or fractal processes f|5], |46|. W. G. Glockle and T. F. 
Nonnenmacher [45[ have formulated fractional differential equations for some relaxation processes 
which are essentially fractal time [HH] processes. 



Fractional Brownian Motion 



Mandelbrot and Van Ness [g4] have used fractional integrals to formulate fractal processes such 
as fractional Brownian motion. The fractional Brownian motion is described by the probability 
distribution B H {t) defined by, 



B H (t) - B H (0) = r(ff 1 + K{t - t')dB{t'), < H < 1, 



(10) 



K(t - 1') 



(11) 



where B(t) is an ordinary Gaussian random process with average zero and unit variance and K(t—t') 

f (t-t') 11 - 1 / 2 0<t'<t 
\ {(t - t') H ~ 1/2 - {-t) H - 1 ' 2 } t'<0. 

The use of the fractional integral kernel is apparent. Recently Sebastian j|7| has given a path 
integral representation for fractional Brownian motion whose measure has fractional derivatives of 
paths in it. 
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Fractional equations for a class of Levy type probability densities 



Levy flights [[48| and Levy walks [|9|, Q have found applications in various branches of physics [[51 
for example in fluid dynamics ]52 and polymers 



In T. F. Nonnenmacher considered a 
class of normalized one sided Levy type probability densities, which provides the length of a jump 
of a random walker, given by 



/(*) 



It is clear that for large x 



r(//) 



/(*) 



X 



exp(— a/x) a > 0, x > 0. 



~ x 



/j > 0, x > 0. 



(12) 



(13) 



where /x is the Levy index. In |^0| it was shown that these Levy-type probability densities satisfy 
a fractional integral equation 



x q t(x) = a H — : 

or equivalently a fractional differential equation given by 

d n -^f\x) d n 



dx n -^ 



dx r - 



x 2 »f(x)). 



(14) 



(15) 



Here, n=lif0</i<l and n = 2ifl<//<2, etc. An interesting observation of this work is 
that the Levy index is related to the order of the fractional integral or differential operator. 

Fractional equations for physical processes 

Modifications of equations governing physical processes such as diffusion equation, wave equation 



and Fokker-Plank equation have been suggested JIT], [12], [54], |55|, |56[ which incorporate fractional 
derivatives with respect to time. In refs. EI, I2"| a fractional diffusion equation has been proposed 
for the diffusion on fractals. Asymptotic solution of this equation coincides with the result obtained 
numerically. Fogedby |57] has considered a generalization of the Fokker-Plank equation involving 
addition of a fractional gradient operator (defined as the Fourier transform of — k^) to the usual 
Fokker-Plank equation and performed a renormalization group analysis. 



Transport in chaotic Hamiltonian systems 

Recently Zaslavasky [^H| showed that the Hamiltonian chaotic dynamics of particles can be described 
by a fractional generalization of the Fokker-Plank-Kolmogorov equation which is defined by two 
fractional critical exponents (oc,/3) responsible for the space and time derivatives of the distribution 
function correspondingly. With certain assumptions, the following equation was derived for the 
transition probabilities P(x,t): 



dPpjx, t) 



d a 
d(—x) 



-{A{x)P{x,t)) + - 



d 



2a 



2d(-x) 



2a 



(B(x)P{x,t)) 



(16) 
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where 

A(r\ = JJtl , 

At-*0 (At) 

and 



,• A'(x;At) , i(TN 
= hm o y (17) 



B(») = lim (18) 
v 7 At^o (At)^ 

with A' and being a th and 2a* 71 moments respectively. The exponents a and /3 have been related 
to anomalous transport exponent. 

We remark that most of the applications reviewed in this section deal with asymptotic scaling 
only. 



3 Weyl fractional differentiability and Holder classes 

The purpose of this section is to review the classical work of Stein, Zygmund and others |59, B^, pTl 
62| which uses Weyl fractional calculus in the analysis of irregular functions. Let f(x) be defined 
in a closed interval I, and let 

u(8) =uj{5: f) =sup\f{x) -f{y)\ for x,y e I \x - y\ < 6 (19) 

The function u(5) is called the modulus of continuity of /. If for some a > we have w(5) < CS a 
with C independent of 5, then / is said to satisfy Holder (In old literature it is called Lipschitz) 
condition of order a. These functions define a class of function A a . Notice that with this definition 
of the modulus of continuity only the case < a < 1 is interesting, since if a > 1, f'{x) is zero 
everywhere and the function is constant. Hence this definition has been generalized (see ||28|| ) to 
one where the case a > 1 is also nontrivial. The idea is to subtract an appropriate polynomial from 
a function / at y. 

Cu(8) = uj(S : f) = sup \f(x) — P(x — y)\ for x,y E I \x — y\ < 8 (20) 

where P is the only polynomial of smallest degree which gives the smallest order of magnitude for 
uj. Recently, this definition has been widely used [12, |H| to characterize the velocity field of a 
turbulent fluid. 

Following Welland [^2| we introduce 

Definition 1 f is said to have an ath derivative, where < k < a < k + 1 with integer k, if D~@f , 
(3 = k + 1 — a, has a k + 1 Peano derivatives at x i.e. there exists a polynomial P xo (t) of degree 
< k + 1 s. t. 

(D-Pf)(x + t)-P xo (t)=o(\t\ k+1 ) t-0. (21) 

Further if 

{- [ P \(D-Pf)(x + t)-P X0 (t)\ P dt} 1/p = o(p k+1 ) (l<p<oo) (22) 

f is said to have an ath derivative in the LP sense. The is in the Weyl sense. 
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It is clear that this definition of fractional differentiability is not local. Particularly the behavior of 
function at — oo also plays a crucial role. The main classical results can be stated using this notion 
of differentiability and involve the classes and N% which are given by the following definitons. 

Definition 2 If there exists a polynomial Q Xo {t) of degree < k s. t. f(xo + t) — Q Xo (f) = 0(\t\ a ) 
as t — > then f is said to satisfy the condition A a and if 

{- f \f(x + t)-Q XQ (t)\ P dt} 1/p = O(p a ), p^O (l<p<oo) (23) 
P J -p 

f is said to satisfy the condition . 



Definition S f is said to satisfy the condition if for some p > 

- r l^° + f >-W dt < oo. (24) 

We are now in a position to state the classical results in the form of theorems 1 to 4. Next two 
theorems pH], give the condition under which the fractional derivative of a function exists. 

Theorem 1 Suppose that f satisfies the condition A a at every point of a set E of positive measure. 
Then D a f(x) exists almost everywhere in E if and only if f satisfies condition N 2 almost everywhere 
in E. 



Theorem 2 The necessary and sufficient condition that f satisfies the condiiton N a almost ev- 
erywhere in a set E is that f satisfies the condition and D a f exisits in the L 2 sense, almost 
everywhere in this set. 

The following results in tell us how the class of a function changes when an operation of 
fractional differentiation is performed. 

Theorem 3 Let < a < 1, f3 > and suppose that f 6 A a . Then D" l3 f 6 A a+I 3 if a + (3 < 1 



Theorem 4 Let < 7 < a < I,. Then D"*f e A Q _ 7 if f e A a 

Though they have their own value, these results are not really adequate to obtain information 
regarding irregular behaviour of functions and Holder exponents. We observe that the Weyl defini- 
tion involves highly nonlocal information and hence is somewhat unsuitable for treatment of local 
scaling behaviour. In the next section we introduce a more appropriate definition. 



4 Local Fractional Differentiability 



In our previous work [RB we introduced the notion of local fractional derivative and demonstrated 
its use in the study of local scaling behaviour. We now briefly explain this notion and use it in 
subsequent sections. 
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4.1 Local fractional derivative and critical order 



Recall the observations made in section [| viz, (1) nonlocal information contained in fractional 
derivative and (2) the fractional derivative of a constant is not zero. The appropriate new notion 
of differentiability must bypass the hindrance due to these two properties. These difficulties can be 
remedied by introducing 

Definition 4 If, for a function f : [0, 1] — > M, the limit 

p^^fflJ!) (25) 
J yyi d(x -y)9 K J 

exists and is finite, then we say that the local fractional derivative (LFD) of order q (0 < q < 1), at 
x = y, exists. 

In the above definition the lower limit y is treated as a constant. The subtraction of f(y) corrects 
for the fact that the fractional derivative of a constant is not zero. Whereas the limit as x — > y 
is taken to remove the nonlocal content. Advantage of defining local fractional derivative in this 
manner lies in its local nature and hence allowing the study of pointwise behaviour of functions. 



This will be seen more clearly in section [4.2| after the development of Taylor series. 

Definition 5 We define critical order a, at y, as 

a(y) = Sup{q | all LFDs of order less than q exist at y}. 

Though these definitions are interesting only when the critical order is less than one, for the same 
reason as that for the first definition (section |3|) of modulus of continuity we extend them for all 
values of a > 0. 

Definition 6 //, for a function f : [0, 1] — > M, the limit 

mM ^ ,(m -^r^-^ pa, 

JKyj *-*v [d(x-y)}i y ' 

exists and is finite, where N is the largest integer for which N th derivative of f(x) at y exists and is 
finite, then we say that the local fractional derivative (LFD) of order q (N < q < N + 1), at x = y, 
exists. 

We consider this as the generalization of the local derivative for order greater than one. The 
definition of the critical order remains the same since, for q < 1, ([25|) and fl26|) agree. This definition 
extends the applicability of LFD to functions where the first derivative exists but are still irregular 
due to the nonexistence of some higher order derivative. As an example we note that the critical 
order of f(x) = a + bx + c|x| 7 at origin, according to definitions 5 and 6 is 7. 
Remark: l)It is interesting to note that the same definition of LFD can be used for negative values 
of the critical order between -1 and 0. For this range of critical orders N = — 1 and the sum in 
equation (|26| ) is empty. As a result the expression for LFD becomes 

E>if(y) = lim r f f{ - % \, (27) 
J yyj *-*y [d(x - y)Y y ' 
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An equivalence between the critical order and the Holder exponent, for positive values of critical 
order, will be proved in section || Here we would like to point out that the negative Holder exponents 
do arise in real physical situation of turbulant velocity field (see [ |63| , |64[ and references therein). 
2) Another way of generalizing the LFD to the values of critical order beyond 1 would have been to 
write it as 

PfM , h jtfM! (28) 

Jyyj *-v [d{x-y)\i v ; 

But the existence of N th derivative of / at x may not be guaranteed in general. Such a situation 
may arise in the case of multifractal functions to be treated in section || 

4.2 LFD and generalized Taylor series 

In order to see the information contained in the LFD we consider fractional Taylor's series with a 
remainder term for a real function /. Let 

F(y, x — y;q, N) = ^—^ (29) 

and 

N f( n )(y) 

F N (x, y) = f(x) - £ ^-^(x - y) n ). (30) 

n=0 1 V 1 "r i ) 

It is clear that 

E>*f(y) = F{y,0;q,N) (31) 
Now, for N < q < N + 1, it can be shown that 

r(g+l)Jo dt v y ; \ j 



i.e. 



where R±(x, y) is a remainder given by 

fl 1 (x,y)= , 1 -fM^-y-^ (34) 

We note that the local fractional derivative as defined above (not just fractional derivative), along 
with the first N derivatives, provides an approximation of f(x) in the vicinity of y. We further note 
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that the terms on the RHS of eqn(^) are nontrivial and finite only in the case q — a, the critical 
order. In ref.|)5j a fractional Taylor series was constructed by Osier using usual (not local in the 
present sense) fractional derivatives. His results are, however, applicable to analytic functions and 
cannot be used for non-smooth functions directly. Further Osier's formulation involves terms with 
negative orders also and hence is not suitable for approximating schemes. 
When < q < 1 we get as a special case 

E> q f(y) 

/(g) = f(v) + p/ + 1 x i x - v) q + Remainder (35) 

provided the RHS exists. One may note in equation (|35|) that when q is set equal to one in the 
above approximation one gets the equation of the tangent. It may be recalled that all the curves 
passing through a point y and having same tangent form an equivalence class (which is modelled 
by a linear behavior). Analogously all the functions (curves) with the same critical order a and the 
same ID a will form an equivalence class modeled by power law x a This is how one may generalize 
the geometric interpretation of derivatives in terms of tangents. This observation is useful when 
one wants to approximate an irregular function by a piecewise smooth (scaling) function. 

5 Fractional Differentiability of Weierstrass Function 

Consider a form of Weierstrass function, viz, 

oo 

W\(t) = \ {s ~ 2)k sin\ k t, A > 1, 1< s < 2, t real. (36) 

k=l 

Note that Wa(0) = 0. The box dimension of the graph of this function is s. The Hausdorff dimension 



of its graph is still unknown. The best known bounds are given by Mauldin and Williams [SB| where 
they have shown that there is a constant c such that 

s < dimngraphf < s. 

log A 

Recently it was shown |67| that if a random phase is added to each sine term in (|36|) then the 



Hausdorff dimension of the graph of the resulting function is s. These kind of functions have been 



studied in detail in PL EH, |69L [70 



We note that there are dynamical systems with graphs of such functions as invariant sets. For 
example, let g : M — >• M be differentiable, and let h : M 2 — > M 2 be given by 

h{x,t) = {Xt,X 2 - s {x-g)t)). (37) 

Then the graph of / can easily be seen to be a repeller of for h, where / is the function given by 

oo 

f(t)^J2X^ k g(X k t). (38) 

k=l 

In the following two subsections we prove lower and upper bound on critical order of the Weier- 
strass function. 
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5.1 Lower bound on critical order 

To check the fractional differentiability at any point, say r, we use t' = t — r and W x (t',r) = 
W x (t' + r) - W x (t) so that W x (0, r) = 0. We have 

oo oo 

W x (t f , r) = ^ {s ~ 2)k sin\ k (t' + r)-J2 \ {s ~ 2)h sin\ k T 

k=l k=l 

oo 

= ^ {s ~ 2)k (cos\ k Tsin\ k t' + sin\ k T(cos\ k t' - 1)) (39) 

k=l 

Now we take fractional derivative of this with respect to t' . 

W ti { d(XH>)<* d{\H'Y ) l4Uj 

From equations (|6]), (|7|) and second mean value theorem it follows that the fractional derivatives 
inside the above sum is bounded uniformly for all values of X k t. This implies that the series on the 
right hand side will converge uniformly for q < 2 — s, justifying our action of taking the fractional 
derivative operator inside the sum. 

Also as t' — > for any k the fractional derivatives in the summation of equation ( 40|) goes to 
zero. Therefore it is easy to see from this that 

d q W\(f r) 

E> q W x (r) = lim , ' ' =0 for q < 2 - s. (41) 

This shows that q th local derivative of the Weierstrass function exists and is continuous, at any 
point, for q < 2 — s. 



5.2 Upper bound on critical order 



For q > 2 — s, right hand side of the equation (J40 ) seems to diverge. We now prove that the LFD 



of order q > 2 — s in fact does not exist. To do this we write the Weierstrass function as follows. 

N 

W x (t) = J2 ^ {s ' 2)k sin\ k t + X^ N W x (X N t). (42) 

k=l 

We now write 

d Q (W X (t) - W X (t')) A Mf^f\k+\ _ „„W\k^ 



» (s „ 2+q)k d*(stn(X k t) - szn(X k t>)) 

2-^ A(\k+\q 



d(t - t'Y ^ d(XH)i 



iy( s-2 + «)N d q (W X (X N t)-W x (X N t')) 
+ d(X N (t-t f ))i K } 



We choose N such that A~ (Ar+1) <\t- t'\ < X~ N . Now since \sin{X k t) - sin{X k t')\ < X k \t - 1'\ 

\t-t'\'- q (44) 



d*\sin{\H) - sin{\H')\ r >*(!_,), 
d[X k (t-t')}i ~ 
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This implies that the absolute value of the first term in equation ( f4H| ) is bounded from above by 
C\( s ~ 2+q ^ N . Now since the first derivative of the Weierstrass function does not exist at any point 
there exists a sequence of points t n approaching t' such that \W x (X N t n ) — W x (X N t')\ > c\ N \t n — 
Therefore there exists a sequence t n such that 

d[X N (t n -t')]i - A ltn 1 [ b) 

and this is valid for every c for large enough n. This implies that 

x (s-2 + a )N d 9 \Wx(\ N tn)-W X (\ N t')\ (s -2 +g) N (m 

d(\»{t n -t'))<l " W 

MW x (t n )-W x (t>)) _ (s _ 2+a)N d«(W x (\ N t n ) - w x (\ N t>)) 

1 d{t n -t'Y ' d{\ N (t n -v)y 1 

< } (s-2+ q )k dq \-^n{^n) - sin(\H')\ ^ 



fc=i 



d(\*(t n - t')Y 



Therefore we get 



This implies that 



^(W X (t n ) - W X (t')) _ cX(s -2 +q )N l < CX (s-2 +q)k 



1 d(t n - vy 



(48) 



(49) 



for C > 0. From this it is clear that the LFD of order greater that 2 — s does not exist. This 
concludes the proof. 

Summarizing, therefore, the critical order of the Weierstrass function is 2 — s at all points. It 
may be noticed that this proof is valid for any A > 1. This generalizes a similar result of @, [36| 
which is valid only for sufficiently large A. Thus there is a direct connection between dimension and 
the differentiability properties for W x (t). As seen below, this observation is not restricted to W x (t) 
but is quite general. 



5.3 Levy index of a Levy flights and critical order 



Schlesinger et al |7T] have considered a Levy flight on a one dimensional periodic lattice where a 
particle jumps from one lattice site to other with the probability given by 



uj — 1 00 

p{x) = ^— u~ j [6{x, +V) + 6{x, -b>)] 



2u 



(50) 



j=Q 



where x is magnitude of the jump, b is a lattice spacing and b > u > 1. 5(x, y) is a Kronecker delta. 
The characteristic function for P(x) is given by 



UJ — 1 00 

P(k) = ^ r -J2uj- j cos(b j k). 



2uj 



(51) 



3=0 
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which is nothing but the Weierstrass function. For this distribution the Levy index is logcu/logfr, 
which can be identified as critical order of P(k). 

More generally for the Levy distribution with index \x the characteristic function is given by 

P(k) = Aexpclfc^. (52) 

Critical order of this function at k = also turns out to be same as \i. Thus the Levy index can be 
identified as the critical order of the characteristic function at k — 0. 



6 Equivalence between critical order and the Holder expo- 
nent 



We recall the definition of Holder exponent h(y) of a function / at y [12], [13]] as h such that there 



exists a polynomial P n (x) of order n, where n is the largest integer smaller than h, that satisfies 

\f( x )-p n ( x -y)\ = 0(\x-y\ h ), (53) 



for x in the neighbourhood of y. When h is restricted between and 1 the equation ([53]) takes the 
form 

\f(x)-f(y)\ = 0(\x-y\ h ), (54) 

The Holder exponent characterizes the behaviour of the function around a given point. One can 
study the function at every point and find its pointwise Holder exponents. According to equation 
(|53l) an analytic function has h — oo at every point. 

As can be seen clearly that the definition of the Holder exponent is not algorithmic and hence 



methods need to be developed for its determination. In |36[ the following result was proved in 
a slightly different form, which establishes an equivalence between the Holder exponent and the 
critical order. 

Theorem 5 The continuous function f(x) has Holder exponent a, < a < 1, at y iff E) q f(y) = 
for all q < a and ID q f(y) does not exist for 1 > q > a. 

This result can be understood from the simple example of fractional derivative of a power x p , given 
in equation (|5[). In this equation if < p < 1 and < q < 1 then the RHS exists. But if we take a 
limit x — > then we get 



d q x p 



if q < p 



lim — — = < const if p — q (55) 
x^o dxi ,f 

oo otherwise 

The LHS in the above equation is nothing but the LFD of x p . This shows that the critical order 
gives the exponent p. This is expected to happen for any function and at any point with a local 
power law behaviour. 

The following corrollary follows from the theorem |5| and a well known result giving relation 
between Holder exponent and box dimension of a graph of a fractal function pi] . 
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Corollary 1 If the critical order of a function f(x) at every point x is a then dirrisf = 2 — a 
where dim^f is a box dimension of the graph of the function f . 

With a slight modification in the proof of theorem |5] a general result giving equivalence between 
the Holder exponent and the critical order using the general definition of LFD follows. The functions 
F and F are defined in equations and ( PUD respectively. 

Theorem 6 Let f : [0, 1] — ► M be a continuous function. 

a) If ID q f(y) = where N < q < N + 1, for some y, then h(y) > q. 

b) If there exists a sequence x n — > y as n — > oo such that 

lirnF(y,x n - y;q,N) = ±oo , 

for some y, then h(y) < q. 

Similarly a following converse of the above theorem can also be proved. 

Theorem 7 Let f : [0, 1] — > M be a continuous function. 

a) Suppose 

\F N (x,y)\ < c\x - y\ a , 

where c>0,N<a<N + l and \x — y\ < 5 for some 5 > 0. Then E) q f(y) = for any q < a for 

z/e (0,1) 

b) Suppose that for y £ (0, 1) and for each 5 > there exists x such that \x — y\ < 8 and 

\F N (x,y)\>c5 a , 

where c > 0, 5 < 5o for some 6q > and < a < 1. Then there exists a sequence x n — > y as 
n — » oo such that 

lim F(y, x n — y; q, N) = ±oo for q > a 

n — >oo 

These two theorems give an equivalence between Holder exponent and the critical order of 
fractional differentiablity. Their proofs are similar to that of theorem |5]. 

7 Isolated masked singularities 

The purpose of this section is to demonstrate the use of LFD to detect masked singularities. We will 
consider only isolated singularities. We choose the simplest example f(x) = J2n=o o. n x n + ax a , N < 
a < N + 1, a;>0. Critical order at x = gives the order of singularity at that point whereas the 
value of the LFD ID q=a f(0), viz aT(ct + 1), gives strength of the singularity. 

Using LFD we can detect a weaker singularity masked by a stronger singularity. As demonstrated 
below, we can estimate and subtract the contribution due to stronger singularity from the function 
and find out the critical order of the remaining function. Consider, for example, the function 

N M 

f( x ) = J2 a n x n + ax a + ]T a n x n + bx^ ', (56) 

n=0 n=N+l 
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where N<a<N + l<M<f3<M + l and x > 0. LFD of this function at x = of the order a 
is D a f(0) = aT(a + 1). Using this estimate of stronger singularity we now write 

rt \ ft \ / (w) (°) n *> a f(0) a 

G(x] a) = f(x) - 2^ x - — — —x . 

^ Q T{n + l) r(a+l) 

The critical order of this function, at x = 0, is (3 which is a masked singularity. Notice that the 
estimation of the weaker singularity was possible in the above calculation just because the LFD 
gave the coefficient of x a /T(a + 1). This suggests that using LFD, one should be able to extract 
secondary singularity spectrum masked by the primary singularity spectrum of strong singularities. 
Hence one can gain more insight into the processes giving rise to irregular behavior. 

Comparison of two methods of studying pointwise behavior of functions, one using wavelets 
and the other using LFD, shows that characterisation of Holder classes of functions using LFD is 
direct and involves fewer assumptions. Characterisation of Holder class of functions with oscillating 
singularity, e.g. f(x) = x a sin(l/x lS ) (x > 0, < a < 1 and (3 > 0), using wavelets needs two 
exponents ||27|| . Using LFD, owing to theorem I and II critical order directly gives the Holder 



exponent for such a function. 

It has been shown in the context of wavelet transforms that one can detect singularities masked 
by regular polynomial behavior [l2j by choosing a appropriate analysing wavelet. (Wavelets with 
first n, for some suitable n, moments vanishing are considered appropriate). If one has to extend 
the wavelet method to unmask weaker singularities, one would then require analysing wavelets with 
fractional moments vanishing. Notice that one may require this condition along with the condition 
on first n moments. Further the class of functions to be analysed is in general restricted in these 
analyses. These restrictions essentially arise from the asymptotic properties of the wavelets used. 
On the other hand, these restrictions are not relevant while using LFD. 



8 Multifractal function 

We saw in section |5| that the Weierstrass function is a fractal function, i.e., it has the same Holder 
exponent at every point. But there are multifractal functions which have different Holder exponents 
at different points. These functions can be used to model various intermitent signals arising in 
physical applications. Since the critical order gives the local and pointwise behavior of the function, 
conclusions of the theorem [|, |6] and |7| will carry over even to the case of multifractal functions where 
we have different Holder exponents at different points. Selfsimilar multifractal functions have been 
constructed by Jaffard [|13[]. We give one specific example of such a function. This function is a 



solution F of the functional equation 

d 

F(x)=Y / ^F(Sr\x))+g(x), x real. (57) 

i=l 

where Si's are the affine transormations of the kind Si(x) = + (with < 1 and fo/s real) and 
Aj's are some real numbers and g is any sufficiently smooth function (it is assumed that g and its 
derivatives have fast decay). For the sake of illustration we choose H\ = [i2 = 1/3, b\ = 0, 62 = 2/3, 
Ai = 3-", A 2 = 3-0 (0 < a < (3 < 1) and 

g(x) = sin(2irx) if x G [0, 1] 
= otherwise. 
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Such functions are studied in detail in [T3|] using wavelet transforms where it was shown that the 
above functional equation (with the parameters we have chosen) has a unique solution F. Further 
at any point F either has Holder exponents ranging from a to (3 or is smooth. A sequence of points 
5^(0), S'i 2 S'j 1 (0), ■ ■ • , Si n - ■ -Si^O), ■ ■ •, where ik takes values 1 or 2, tends to a point in [0, 1] (in 
fact to a point of a triadic Cantor set) and for the values of /^s we have chosen this correspondence 
between sequences and limits is one to one. The solution of the above functional equation is given 
by m 



-1 



V(*))- 



(58) 



n=Q 



Note that with the above choice of parameters the inner sum in (|58|) reduces to a single term 



Jaffard [|T^] has shown that the local Holder exponent at y is 



h(y) = liminf , 



\ x in{y) l 



(59) 



where {ii(y) • • • i n (y)} is & sequence of intergers appearing in the sum in equation(|58"l) at a point 
y. It is clear that h m i n = a and h max = [3. The function F at the points of a triadic cantor set 
have h(x) G [oc, /3] and at other points it is smooth ( where F is as smooth as g). Benzi et. al. JF2 
have constructed multifractal functions which are random in nature unlike the above nonrandom 
functions. For still another approach also see [[73 



Several well-known 'pathological' functions have been reanalyzed in 3U, 3T, [£J and found 
to have multifractal nature. Here we consider one example of classical multifractal function. 



R{ x ) — E — 2 s ^ n ( 7rn2x 

n=l 71 



(60) 



This function was proposed by Riemann. It turns out that the regularity of this function varies 
strongly from point to point. Hardy and Littlewood |74| proved that R(x) is not differentiable at xq 
if xq is irrational or if xq can not be written as 2p + l/2q+ 1 (p,q G N). In fact they showed that the 
Holder exponent at these points < 3/4. Gerver |75| proved the differentiability of R(x) at points of 
the form 2p + l/2q + 1 (p, q G N). At these points the Holder exponent is 3/2. This function has 
also been studied in [|76l [77|. Jaffard |23] has recently shown that the dimension spectrum of the 



MM- Jaffard 
Riemann function is given as below. 

d(a>) = 



4a - 2 if a G [f , f] 
if a = | 

— oo otherwise 



(61) 



where d(a) gives the Hausdorff dimension of the set when the Holder exponent is a. 

LFD forms one method of studying pointwise behaviour of such multifractal functions alongwith 
other methods and may cosiderably reduce the analysis involved. This fact was demonstrated in 
on a specific example of self-similar multifractal function given by equation (j58|). 
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9 Concluding remarks 



First we reviewed the applications of various fractional differential equations in different physical 
situations. It was noted that most of these applications dealt with the asymptotic scaling. Further 
we reviewed the classical results within the framework of Weyl fractional calculus. However these 
results were found to be inadequate for the study of pointwise behavior of fractal and multifracal 
functions. The notion of LFD as developed in was found suitable for this purpose. In terms of 
the LFD it was also possible to write a fractional Taylor series of a function (useful in an analytic 
treatment of approximations). It was pointed out that generalization of the notion of tangents to 
the graph of a function (useful for geometric purposes) is also possible using LFD. 

It was established that the critical order of the Weierstrass function is related to the box di- 
mension of its graph. If the dimension of the graph of such a function is 1 + 7, the critical order 
is 1 — 7. When 7 approaches unity the function becomes increasingly irregular and local fractional 
differentiability is lost accordingly. Thus there is a direct quantitative connection between the di- 
mension of the graph and the fractional differentiability property of the function. This is one of 
the remarkable conclusions of the new approach. An important consequence of this approach is 
that a classification of continuous paths (e.g. fractional Brownian paths) or functions according to 
local fractional differentiablity properties is also a classification according to dimensions (or Holder 
exponents) . 

Also the Levy index of a Levy flight on one dimensional lattice is identified as the critical order 
of the characteristic function of the walk. More generally, the Levy index of a Levy distribution is 
identified as the critical order of its characteristic function at the origin. 

We have argued and demostrated that LFDs are useful for studying isolated singularities and 
singularities masked by the stronger singularity (not just by regular behavior). It was also shown 
that the pointwise behavior of irregular (fractal or multifractal) functions can be studied using the 
methods of this paper. 

We note, however, that the treatment of random irregular functions as well as multivariable 
irregular functions is badly needed. We hope that these problems can be tackled in near future. 
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